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Quasi two-dimensional pattern forming systems with 
spontaneously broken isotropy represent a novel symmetry 
class, that is experimentally accessible in electroconvection of 
homeotropically aligned liquid crystals. We present a weakly 
nonlinear analysis leading to amplitude equations which cou- 
ple the short-wavelength patterning mode with the Goldstone 
mode resulting from the broken isotropy. The new coefficients 
in these equations are calculated from the hydrodynamics. 
Simulations exhibit a new type of spatio-temporal chaos at 
onset. The results are compared with experiments. 
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Ac driven electroconvection (EC) in nematic liquid 
crystal (NLC) layers is one of the richest systems for the 
study of pattern forming phenomena . We consider 
the typical thin-layer geometry with the (slightly con- 
ducting) NLC sandwiched between glass plates coated 
with transparent electrodes. An ac voltage U is applied 
across the electrodes. By appropriate surface treatment 
of the glass plates the director n (the preferred orienta- 
tion of the NLC molecules) can be fixed at the bound- 
aries. In particular the case of planar alignment (n in the 
plane of the layer) has been studied intensely. At onset 
one may then have normal rolls, where the wavevector is 
parallel to the undistorted director, or oblique ones. 

The case of homeotropic surface anchoring (n aligned 
perpendicularly to the boundaries) where the system is 
isotropic in the plane of the layer (= x-y plane) offers 
novel possibilities. Then, in the traditional EC materials 
with negative dielectric anisotropy, the voltage applied 
across the layer will first turn the director away from 
the layer normal (bend Freedericksz transition) leading 
to a quasi-planar director configuration (see e.g. ]^)). The 
spontaneously chosen direction of the bend (i.e. direction 
of the projection of n on the x-y plane) will be denoted 
by c — c(r, t) ( |c| = 1 ). Further increase of the voltage 
will eventually generate EC in close analogy to the pla- 
nar case Q; in fact nucleation to normal, oblique and 
traveling rolls has been observed S. The notable dif- 
ference to the planar case is that the preferred axis (the 
in-plane director c) is degenerate and not fixed exter- 
nally (neglecting unavoidable inhomogeneities and in the 
absence of a planar magnetic field). Then oblique rolls 
are not expected to lead to a stable pattern because they 
will in general exert a torque on c which cannot be com- 
pensated. Even normal rolls, where a torque is absent 



because of symmetry, can be unstable because transverse 
modulations can be enhanced by the torque. 

Here we will investigate the situation by setting up 
a novel weakly nonlinear description that incorporates 
the critical convection mode together with the Gold- 
stone mode resulting from the spontaneous breaking of 
the 0(2) symmetry by the bend Freedericksz transition. 
The general form of the equations is derived from symme- 
try considerations. Analyzing the stability of rolls indi- 
cates that one may expect spatio-temporal chaos (STC) 
at onset under very general conditions. In fact this is 
probably the first experimentally accessible example of 
a nonrotating system with a direct transition to STC in 
a stationary (non-Hopf) bifurcation J5]]6). The new co- 
efficients of the equations are then calculated from the 
underlying standard hydrodynamic theory making quan- 
titative comparison with experiments on MBBA (JV-(p- 
methoxybenzylidene)-p-butylaniline) without adjustable 
parameters in principle possible. We end up by showing 
results of simulations. 

First we consider the situation where the orientation of 
c differs only by a small angle tp from an overall direction 
that we choose as the x axis (then c = (costp, simp, 0) ~ 
(X,tp,0)). For simplicity we present the analysis for the 
normal-roll case, but it was done also for oblique rolls and 
some results will be given for that case. Let q c = (g c , 0), 
be the wavevector of the most unstable convection mode 
when ip = 0. Near threshold the physical state can be 
written in the form 



V(r,z,t) = ( V^(z,t)A(r,t)e zcl ^ r + c.c.) 
+ V G (z,t)tp(r,t) + h.o.t., 



(1) 



where V is a formal vector consisting of the three com- 
ponents of the velocity, the director n, and the charge 
density. A = const, and <p = yields the critical mode at 
the threshold = e = (U 2 - U?)/U 2 . The reduced de- 
scription in terms of the modulation amplitudes A and 
tp has to respect the symmetries of the underlying hy- 
drodynamic equations and the Freedericksz state. Thus 
we require invariance of the equations under translation 
(r — > r + r , A — > Ae~ jqc ' r °), rotation (tp — > p + S, r — > 
r+<5zxr, 8 X -> d x -6d v , d y -> d y +Sd x , A -» Ae m ****), 
three dimensional inversion (x — > — x, y — > —y, A—* A*), 
reflection at the x axis (y — > —y, ip — > — tp), and time- 
translation (t — > t + to). We consider the regime near 
onset (e, |A| <C 1) where slow spatial and temporal vari- 
ations (d Xl d y , dt <C 1) can be assumed. However, similar 
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to |t[] no specific scaling relations are assumed at this 
point. 

By retaining only the leading terms allowed by sym- 
metries the equations describing the evolution of A and 
tp can be determined up to prefactors: 



rd t A 



g\A\ 2 +il3 y <p, 



A. 



(2) 



(3) 



%d t ip = K x d 2 y y + K 3 dl<p - x a H 2 cp 

r 

+ j(-iqcA*(d y - iq c <p)A + c.c). 



In (|2|) and in similar equations the derivative operators 
d x , d v operate on A only, while ip tV is the derivative of ip 
in the y direction. All coefficients are real and we will as- 
sume all but (3 y , r to be positive. The coefficients (3 y and 
r describe the novel effects. The last expression in Eq. (||) 
describes the torque exerted by the roll pattern on c. An 
additional term —x a H 2 ip is included in Eq. (||) to allow 
a continuous transition from the rotationally invariant 
case to the case of fixed c, as for planar ly aligned cells. 
Physically H corresponds to a stabilizing magnetic field 
parallel to the x-axis. In the special case f3 y — —£, yy q c 
the angle <p would play the role of a gauge field and Eq. 
(0) could be derived from a potential. If also T > 0, 
both equations can be derived from a Lyapunov func- 
tional. However, in our non-equilibrium system there is 
no reason to assume existence of a potential. For MBBA 
we find r < 0, which implies a "repulsive torque" be- 
tween the direction of the local wavevector and c. Then 
at H = all ordered roll solutions are in fact unstable 
(see below). 

For H = Eqs. (|||) become uniform in e when the 
scaling d x ,d y ,A,tp ~ e 1 / 2 and d t ~ e is adopted. This 
is a major difference to similarly looking coupled ampli- 
tude equations describing mean-flow effects in isotropic 
Rayleigh-Benard convection Since ip is also scaled, 

this is not applicable to Eqs. (|l3|,14|) below because there 
(f is 2n periodic. 

After rescaling (@,ph one arrives at the normal form 



fd t A 



l + dl + {dy- iiff 

\A\ 2 + ij3 v if iy A, 



d t Cp = c>l<p + k z dl<p-e- Y H 2 q> 
+ f (-iA*(d y - iip)A + c.c), 



(4) 



(5) 



where A = (e/g) 1/2 A, x = £ xx e 1/2 x, y = £, yy e 1/2 i 

-1/2 - l(c n \ + — X,. C2 f/(rTC\ h — a l(& 



ip = S^v/tiyyqc), t = Myyt/ieKx), f3 y = (3y/(% y q c ), 



K 3 



K 3 e vv /(Kie xx ), H 2 



XattyH 2 /^ f 



qUlyT/^gKi), f = ^ir/( 7l e 2 a ). The caret has been 
suppressed on the variables x, y, t in Eqs. (T 



The stationary, uniform solutions of these equations 
are A = A e i ^ x+p y\ <p = ip = 2fA 2 P/{H 2 /e + 2T A 2 ), 
Aq = I — (P — Lpo) 2 — Q 2 ■ We now present the main re- 
sults of the stability analysis of these roll solutions. For 
f < and H = rolls are always unstable at onset, 
and one expects a (direct) transition to spatio-temporal 
chaos. For H ^ and y < one has a kind of Eckhaus 
instability associated with a finite angle between the di- 
rector and the local wavevector of the rolls, whereas for 
f3 y > an instability induced by director modulations 
via the /3^-term in (^]) is dominating. 

Perturbations can in general be cast in the following 
form: 

§A = (A + e st+i ( k * x+k y y ') + a~ e s " t ~ i ( k * x+k yy)'} e i (Q x + p y) 
Sip = (^ e st+l(fe ^ +fe ^ ) + c.c). (6) 

The case ipo = = P is particularly simple. Then the 
constant-amplitude solutions are stable if H 2 je + 2T(1 — 
Q 2 ) > 0, H 2 /e + 2f (I - Q 2 )(l + $ y ) > and Q 2 < 1/3 
(Eckhaus criterion). For H = stable rolls require T > 0, 
1 + (3 y > and Q 2 < 1/3. Otherwise one has stability if 
PT/ipo > 0, Pyipo/P > -1 and detM^ > 0, where 



M xx = 1 - (P - ^ ) 2 (1 + 2p /P) - 3Q 2 



M 



xy 

Myy = (1 + V P O /P)(1 - 3(P - P ) 2 - 



M„ 



(P-<p )Q((l-j3 y )<p /P-2), 



(7) 



The direction of (k x ,k y ) cr it at the marginal point 
det Mij = is given by the null space of My . 

At the generalized Eckhaus boundary (P — ipo) 2 + 
Q 2 = 1/3 known for planarly aligned EC Q, detM^ ~ 
— (ipo/P) 2 is already negative. Thus, as expected, the 
additional degree of freedom ip leads to a destabilization 
of the system. 

Expanding for small P and Q we find the uniform so- 



lution to be stable for H 2 > H 2 lit where 



= -2TAI 



,1/2 



PI +3P 2 +h.o.t 



(8) 



if (3 y < and 



#crit 



-2fi 2 



l + Pv + 



P 2 Q 2 + h.o.t. (9) 



if (3 y > 0. Since in the latter case the r.h.s. is dominated 
by the quadratic dependence of Aq on P, Q, small but 
finite P, Q are more stable then P, Q = 0. 

We expect that an experimental verification of the sta- 
bility analysis is possible. Note that the value of P can 
be adjusted by turning the magnetic field with respect to 
the probe after a stable roll pattern has been established. 

The analogous equations for the amplitudes A\, A^ 
of oblique rolls (zig and zag) have additional terms and 
the expression for the torque on c acquires the form 
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Tq c p c (\Ai\ 2 - \A 2 \ 2 )/2. If the crossed- roll solution | Ai| = 
\A 2 \ is preferred over simple rolls (e.g. A\ ^ 0, A 2 = 0) 
the situation is presumably analogous to that of normal 
rolls (no torque on the director for the unmodulated so- 
lution) . In the opposite case (as for MBB A) the resulting 
equations have no nontrivial time-independent roll solu- 
tion. In that case e cannot be scaled out completely and 
one gets different scaling ranges. 

Next we describe the calculation of the coupling coef- 
ficients P y and r from the standard hydrodynamic equa- 
tions for NLCs 0,0. The other coefficients are ob- 
tained following standard procedures (see e.g. 
As usual the set of equations is decomposed into a sum of 
operators operating on the state vectors V = V(r, z, t). 
Each operator is linear in each of its arguments: 



= O(V) = -BdtV + LV + N 2 (V\V) 



(10) 



The operator O(V) depends on e in two ways: by the ex- 
plicit action of the electric field and the field dependence 
of the Freedericksz state, which is the basic state in this 
formalism. 

Two components can be distinguished in the kernel of 
L: the convective mode V^ 1 ) at wavevectors q = ±q c 
and the Goldstone mode Vg at q = Q normalized 
such that it corresponds to Sn y — dc/ dip\ v= $. Near 
onset these two modes govern the dynamics described 
in Eqs. (H,||). We will also need the linear eigenmodes 
V (1) (p), V G (p) of © into which and V G de- 

velop for transverse modulations, i.e. at q = ±(q c ,p) 
and ±(0,p) respectively. The coefficients are obtained in 
some analogy to the treatment of mean-flow effects |?]-||] 
and of the concentration field in binary fluids . 

Let (Vi,V 2 ) be a scalar product and (p) and 

Vg + be the adjoint eigenvectors of L corresponding to 
VW(p) and V G . The coefficient (3 y in Eq. (§) turns out 
to be 

/j y =-(v( 1 ) + ,(a e L)v( 1 ))- 1 ^vw + b), (ii) 

(N 2 (V« I V g (p) ) + N 2 ( V G (p) | V« ) 



p=0 
e=0 



The coefficient of the hydrodynamic angular momentum 
in (i) is 



2 



dp 



,N 2 (V(p)\V(p))) 



p=0 



(12) 



where C is chosen such, that 71 = 1. 

In order to explore the attractor in the roll-unstable 
case and to compare the long-time dynamics with exper- 
iments the assumption that if be small must be relaxed 
(at least at H = 0). One in fact only needs the condi- 
tion, that the local wavevector differs little from a (local) 
critical wavevector q c c. Formally this can be achieved 
by reintroducing the fast variation (A = Ae !q <=' r ) and 



rewriting derivatives and the magnetic torque in rotation- 
invariant form (we write 6j_ = z x c): 



rd t A 



e + C x {c ■ d r - iq c ) 2 + ^Jc ± ■ d r ) 



(13) 



g\A\ 2 + i/3 y c ± ■ \Jip 



A. 



71 d t ip = K 3 c ± • V 2 c + [K x - X 3 )c ± • V(V • c) (14) 

' ~*(c± ■ d r )A+c.c). 



r 

+ X a{c • H)(c ± • H) + -(-iq c A*(c A 



These equations have been used for numerical simula- 
tions. 

In Fig. [j] we compare snapshots of ip = A + c.c. from 
our simulations (left side) with experiments (right side) 
from HH ( £ = 0.02, H = 0). Material parameters as in 
were used to calculate the coefficients of ( [l3|]l4"| ). For 
normal rolls (upper pair), at the experimental driving 
frequency ujtq = 1.75 (to = eoe±/ a± is the charge relax- 
ation time) the calculated coefficients are: r = 0.0292, 
q c = 2.71, e xx = 0.434, e yy = 0.225, (3 y = -0.60, 
K x = 1.01, K 3 = 1.28, g = 4.37, T = -12.8, X a = 1- 
(We measure length in units of d/n with d the thick- 
ness of the fluid layer, time in director relaxation times 
Td = 7id 2 /fcn7r 2 and magnetic fields in Freedericksz fields 
H 2 = kuTT 2 /xad 2 -) For oblique rolls (at lut q = 0.73, 
lower pair) we proceeded analogously. In order to obtain 
a good fit to the wavelength the thickness d — 28.5 /j,m 
was used instead of the nominal value d = 20 /im given 
in |l4|] (d was not measured directly) . This only rescales 
length. 

The patterns in experiments and simulations look very 
similar. Though in the normal roll regime some de- 
fects appear, the rolls are locally aligned along some 
main direction. In the simulations we find the pattern 
to change on a time scale t c /e = 0.8 (t c is defined by 
(ip(r,t- to)ip(r,t )) Ttto ~ exp(-|t|/i c )) while the exper- 
imental state appeared to be time independent. Presum- 
ably some pinning of the pattern in the experiments is 
responsible for this disagreement. The oblique roll regime 
is dominated by a superposition of zig and zag. Again 
the preferred axis changes only over large distances. A 
persistent time dependence was observed in simulations 
(to/e = 0.5) and experiments. Unfortunately the exper- 
iments are still preliminary so far. A systematic experi- 
mental investigation is under way p5| . 

Equations (|^^) represent normal-form equations for 
quasi-2D pattern-forming systems with a novel kind of 
symmetry, and should thus be of general interest. Other 
realizations might be found in convection instabilities in 
smectic-C liquid crystals, where one has a c-director from 
the beginning on, and in Rayleigh-Benard convection of 
homeotropically aligned NLCs with an additional electric 
field. In this case the fields that drive the Freedericksz 
transition and the convection instability, respectively, can 
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be varied independently, which should allow to access a 
large parameter range of Eqs. (0jq). 

For H = and T < Eqs. (fb||lj) describe STC at 
onset. There is some analogy to a recent investigation 
of a Id model for seismic waves, where a spontaneous 
symmetry breaking is also important p6| . Other exam- 
ples of STC at onset are the Kiippers-Lortz instability 
(in rotating Rayleigh-Benard convection |l7],[l^]) and sys- 
tems undergoing a Hopf bifurcation. Here one has the 
Benjamin-Feir destabilization mechanism (see e.g. JTg| ] 
and as a recent example |2(|) and the so-called disper- 
sive chaos where a quantitative description by a simple 
Ginzburg-Landau equation should be possible Ell]. The 
origin of chaotic behavior is of course very different in the 
various systems and their detailed comparison appears 
fruitful. Work on the characterization of the complex 
spatio-temporal states is in progress, with particular em- 
phasis on the transition to order under the influence of 
the stabilizing magnetic field. 

For ordinary stationary bifurcations (and Hopf bifurca- 
tions in the Benjamin-Feir stable range) the onset of STC 
is presumably always at finite e (possibly quite small) , at 
least as long as effects from lateral boundaries play no 
role. Examples are the usual case of EC with planar 
alignment p2fl , besides our system in the presence of a 
stabilizing magnet field. 
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[ p3| (right) for normal (top) and oblique (bottom) rolls (see 
text). 
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